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ABSTRACT 

In this letter, we exploit the gravitational potential of the rich cluster A 2218 as a magnifying glass. We 
demonstrate that the magnification due to the cluster allows us to observe distant background galaxies 
at a comparable level of detail to galaxies at z ^ 0.1. Using the GMOS Integral Field Unit on Gemini 
North we observed the spatially-resolved [On] A3727 emission line spectrum for a lensed disk-galaxy 
at z=1.034. Using a detailed model for the cluster mass distribution, we are able to correct for the 
lensing by the cluster and reconstruct the source morphology. We find that the overall magnification 
is a factor of 4.92 ± 0.15, and the rest-frame absolute /-band magnitude is Mj*^** — —22.4 ± 0.2, where 
the error bars include conservative estimates of the uncertainty in the source-plane reconstruction. The 
inclination-corrected circular velocity is 206 ± 18kms~^. The galaxy lies very close to the mean Tully- 
Fisher relation of present-day spirals. Although our results are based on a single object, they demonstrate 
that gravitational lensing can be viably used to make detailed studies of the evolution of the structure 
of distant field galaxies. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: halos — galaxies: 
high-redshift — galaxies: kinematics and dynamics — galaxies: spiral 



1. INTRODUCTION 

The defiection of light-rays by the deep gravitational 
potential of galaxy clusters can be harnessed to greatly 
increase the effective collecting area of astronomical tele- 
scopes. Such "gravitational telescopes" magnify the im- 
ages of background objects allowing us to study faint dis- 
tant galaxies in a level of detail that would simply not be 
possible by conventional means (Small et al. 1996; Franx 
et al. 1997; Teplitz et al. 2000; EUis et al. 2001; Campu- 
sano et al. 2001; Smith et al. 2002). In particular, many of 
the selection biases that require us to only observe bright, 
extended galaxies are removed, allowing us to test models 
for galaxy formation in a much fairer way than otherwise 
possible. 

In this paper, we will concentrate on the dynamics of 
galaxies seen at z 1 by measuring the rotation speeds of 
gas in these systems through the [On] A3727 emission line. 
By observing the distant (z > 1) universe, we sample cos- 
mic history at a time before many present-day stars had 
formed. The "classical" galaxy formation model, where 
the formation of a dark halo preceeds the formation of the 
galaxy's disk (with the disk subsequently growing over an 
extended period of time through the gradual accretion of 
gas from the halo, eg. Eggen et al. 1962; Larsen et al. 1980; 
Sandage 1990), predicts that the stellar masses at a given 
circular velocity will be much lower at z = 1 (by about 
a factor of 2) than at the present-day. In contrast, CDM 
hierarchical models predict that galaxies would follow a 
similar stellar mass versus circular velocity correlation at 
all redshifts. Previous work in this area has, by necessity, 
concentrated on galaxies which are intrinsically bright and 



et al. 2002; Ziegler et al. 2003), with only a few galaxies 
being studied at higher redshift (Vogt et al. 1996; Milvang- 
Jensen et al. 2002; Barden et al. 2003). The results suggest 
an increase in B-band luminosity for a given circular veloc- 
ity (Vogt et al. 2002; Kannappan et al. 2002) - suggesting 
a preference for hierarchical formation models. 

Here we demonstrate the feasibility of using lensing to 
study the evolution of the TuUy-Fisher relation at much 
higher redshift. We use the GMOS IFU on Gemini-North 
to target the gravitationally lensed z = 1.034 arc (object 

# 289 - we use the same numbering system as Pello et 
al. (1992) in the cluster A 2218. We use a cosmology with 
Ho = 75kms~^Mpc~^ and qo = 0.05 to be consistent with 
Vogt et al. (1997). We note that for a cosmology with 
VIm = 0.3, JIa = 0.7, the results are very similar (the 
physical distances are only 6% larger and the luminosities 
12% brighter). 

2. OBSERVATIONS, ANALYSIS & RESULTS 

The galaxies lensed by the A 2218 cluster were first stud- 
ied by Pello et al. (1988), and later a number of arcs and 
arclets were discovered in the cluster core (Kneib et al. 
1996; Ebbels et al. 1998; EUis et al. 2001). Spectroscopy 
of the arcs in the cluster have confirmed a redshift for arc 

# 289 of 1.034 (Pello et al. 1992; Ebbels et al. 1998). 

2.1. HST Observations and Lens Model 

The cluster was first observed by the Hubble Space Tele- 
scope (HST) WFPC2 in September 1994 in the R702 fil- 
ter for a total of 6.4 ks (see Kneib et al. 1996). It was 
subsequently observed again in January 2000 using the 



mostly at z < 0.5 (Vogt et al. 1997; Verheijen 2001; Bohm 

^ Department of Physics, University of Durham, South Road, Durham DHl 3LE, UK - Email: a.m.swinbank' 
^ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS OHA, UK 
3 California Institute of Technology MC 105-24, Pasadena, CA 91125, USA 
* Observatoire Midi-Pyrenees, 14 Avenue E. Belin, 31400 Toulouse, France 

1 



idur.ac.uk 



2 






Fig. 1.— Left: Arc #289 in A2218 generated by combining the HSTWFPC2 5450^606 and /gw drizzled images. (Middle) The [On] A3727 
emission map of the arc measured from our IFU observations. The distribution of [On] emission agrees well with the UV flux seen in the left 
panel (the seeing for the observations was 0.7"). (Right) The contour map of the derived velocity field of the galaxy. The scale is marked in 
aJTcseconds, North is up and East is left. 



B45o,V606 and I8i4 filters for totals of 12.0, 10.0 and 12.0ks 
respectively (discussed in Small et al. 2001)^. 

Arc#289 is a fairly blue disk-galaxy at a = 
16''35'"55M, S = +66°11'51".0 (J2000) lying in a sad- 
dle between the central dominant galaxy of the cluster (# 
301) and a sub-clump of elliptical galaxies 307 & # 
292, Figure 1 in Small et al. 2001). The if5T images of 
the galaxy show a large amount of internal structure with 
the knotted, disk-like morphology resembling a late-type 
galaxy (especially the B-band flux which samples the rest- 
frame UV and is therefore dominated by the star-forming 
Hll regions - Figure 1 - left). We calculate the magnitude 
of the arc in various passbands in Table 1 by using the 
IRAF IMSURFIT package with the sky estimated from a 2"*^ 
order polynomial surface fit to regions around the galaxy. 
We then use SEXTRACTOR (Bertin & Arnaud 1996) to es- 
timate the residual background within the frame. 

Table 1: Photometric properties of #289 



reconstruct) the detailed properties of high redshift galax- 
ies. 

While the giant blue arc # 289 is highly distorted at the 
northern end (Fig. la), extending across the halo of the 
cluster galaxy # 244, the majority of the source lies outside 
the caustic producing a highly magnified, weakly sheared 
image. The southern end of the arc does not suffer the 
same strong shear, and is not so distorted although it has a 
higher surface brightness than the northern end. Exploit- 
ing the fact that gravitational lensing conserves surface 
brightness, we compute the inagnific;atioii in each band 
by using the transformation between the sky-plane and 
source-plane coordinates and compare the flux in the (re- 
constructed) source and (observed) sky frames (Fig. 2 - 
left). We use Smith's (2002) model of A 2218 to compute 
the magnification (n) of # 289, obtaining a mean lumi- 
nosity weighted magnification of /U = 4.92 across the whole 
galaxy image. The spatial variation of surface brightnesses 
may vary between filters due to differing spatial distribu- 
tions of stellar populations (and dust) in the lensed galaxy, 
thus causing the magnification to vary by ±0.15 depending 
on the band. The statistical error bar on fi is derived from 
the family of lens models which adequately reproduce the 
multiply images arcs in the cluster. We estimate the 1-a 
uncertainty by perturbing the parameters of the best fit 

Aperture magnitudes of the arc are in the Vega-based system. lenS model SUCh that = 1. For each model we reCOm- 

inciination is computed in the source (reconstructed) frame pute the magnification of #289. The error Corresponds to 

the largest variation in magnification that we found and 
To correct for the distortion and magnification of the jg converted into a cjonservative estimate of ±0.05 magni- 
galaxy image by the cluster lens we need to employ a lens tudes in the source plane photometry of #289. 
model. A detailed mass model for A 2218 was originally 
developed by Kneib et al. (1996) from HST imaging and 
this has been updated by Ellis et al. (2001) and Smith et 
al. (2003). In particular. Smith (2002, see also Smith et al. 

2003) incorporated all of the strong lensing constraints to ^ ^ ^ , , , t ■ 

J r.u u J- 4. • J 4. 1 • 1 J. 2.2. Ground-based Imaqinq 

produce one or the best constrained strong lensing clusters, " " 

with three spectroscopically confirmed multiple-image sys- Additional constraints on the spectral energy distribu- 

tems. The precision of this lens model makes this cluster tion (SED) of the galaxy comes from the near-infrared 

and ideal gravitational telescope with which to study (and imaging of this field using the INGRID camera (Packham 

^ Based on observations made with the NASA/ESA Hubble Space Telescope which is operated by STSCI for the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Magnitude 


b/a 


inclination 


-B450 


21.66 ± 0.04 


0.37 


68 


^606 


20.86 ± 0.10 


0.41 


66 
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20.53 ± 0.04 


0.40 
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^814 


19.79 ± 0.05 


0.53 
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18.63 ± 0.06 
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et al. 2003) on the WiUiam Herschcl Telescope 2. The re- 
duction and analysis of these data is described by Small 
et al. (2001). With our multi-wavelength photometry (Ta- 
ble 1) we construct an SED of this very blue galaxy whose 
colors indicate a current to past star formation rate con- 
sistent with that seen for late-type spirals at the present 
day. The rest frame /-band luminosity is approximately 
equivalent to the observed i/-band and so is calculated 
from the J and K photometry (Table 1). We co-add the 
J. K images to obtain an aperture which is used to ex- 
tract the J, K magnitudes and interpolate for the relevant 
SED type. We derive a rest-frame /-band magnitude of 
Ml = —24.1 ± 0.1. Correcting for the magnification due 
to the lens, the absolute rest-frame /-band magnitude is 
y^rest ^ -22A ± 0.2. We apply the same technique for 
the rest-frame B-band magnitude, (at 2; ~ 1, / and R- 
bands are the closest match to rest-frame B-band) and 
estimate its uncertainty by computing the magnitude for 
a variety of SED types which are consistent with the ob- 
served optical colours. We compute a corrected rest frame 
S-band magnitude of M|j«^* = -21.1 ± 0.2. Accounting 
for the (1 + zY dimming, the central surface brightness of 
the galaxy {i-ib,a ^ 20.6) confirms #289 as a high-surface 
brightness (HSB) galaxy. 

The inclination is determined by transforming the HST 
images of the galaxy to the source plane. We fit ellipses 
to an isophote of the galaxy image in the source plane us- 
ing the IDL GAUSS2dfit routine. The elliptieity is then 
e = 1 — b/a (where a and b are the major and minor axis 
of the ellipse) and the inclination, i, is cosi = b/a. The 
average axis ratio of the ellipses in the four passbands (Ta- 
ble 1) is 0.43 ± 0.07 which translates to an inclination of 
64 ± 6 degrees assuming an intrinsically circular disk. 

101 I 101 ' ' ' 
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Fig. 2. — The reconstructed image of the arc corrected for lens 
magnification using the mass model of Smith et al. (2003). Left; the 
reconstructed image of the galaxy based on HST imaging. Right: 
The velocity map of the galaxy in the source frame. The dark 
and light regions represent redshift and blueshift respectively and 
the contours map the velocity in the source plane. The luminosity 
weighted magnification of the source is 4.92 ± 0.15, but varies from 
~ 5.6 to ~ 4.9 from the northern to southern end of the arc. The 
scale shows the size of the galaxy in the source frame. Without a 
lens, at z = 1, 1" corresponds to 7.7kpc. 



2.3. GMOS Spectroscopic Imaging 

A2218-arc # 289 was observed with the GMOS-IFU 
on Gemini North on 2002 .June 12'^^ U.T. during Science 
Demonstration time for a total of 5.4 ks in 0.7" seeing 
and photometric conditions. The IFU uses a lensed fiber 
system to reformat the 7" x 5" field into two long slits 
(AUington-Smith et al. 2002). Using an /-band filter in 
conjunction with the R400 grating results in two tiers of 
spectra recording a maximum field of view. The spectral 
resolution of this configuration is A/ A A = 2000. For the 
galaxy at z = 1.034, the emission for the [On] doublet falls 
at a wavelength of 7581 A, in a region of low sky emission. 

The GMOS data reduction pipeline was used to extract 
and wavelength calibrate the spectra of each IFU element. 
The variations in fiber-to-fiber response were removed in 
IDL by using continuum regions either side of the expected 
range of [On] emission. The [On] doublet is clearly re- 
solved in the GMOS spectra (Fig. 3). This is useful, since 
one of the [On] lines can clearly be identified even when 
the other lies immediately under the narrow sky line at 
^ 7580A. The emission line doublet was fitted using a 
minimisation procedure, taking into account the greater 
noise at the position of the sky line. The spectra were 
averaged over a 3 x 3 spatial pixel region, increasing this 
region to 4 x 4 pixels if the signal was too low to give a suf- 
ficiently high improvement over a fit without the line. 
In regions where this averaging process still failed to give 
an adequate x^) no fit was made. With a continuum fit 
we required a minimum of 25 (S/N of 5) to detect the 
line. To compute the error, the velocity parameter of the 
best fit line is varied until the signal drops by a of 9. 
The range in velocity corresponds to a formal 3(7 error. 

We use the velocity field to infer the rotational velocity 
of the galaxy's gas disk. As can be seen from the major 
axis cross section, (Fig. 4 inset) the asymptotic rotation 
speed is Vrot = 186±16kms~-^, (i.e. the total velocity shift 
across the galaxy is 2vrot)- The inclination corrected ve- 
locity of the galaxy is 206 ±17 kms"""^. This corresponds 
to a dynamic mass of 1.4 x IO^^MqHj^ within a radius 
^ lAh^^kpc in this cosmology. 

In addition the shape of the rotation curve can be com- 
pared to the low redshift {z ~ 0.03) sample of spiral galax- 
ies by Courteau (1997). We fit the observed source plane 
1-D rotation curve by convolving the ARCTAN and multi- 
parameter models (Courteau 1997) with 0.7" seeing (trans- 
formed to the source plane). Using a x^ fit the observed 
rotation curve is best described with an arctan function 
with transition radius, rt = 1.7kpc (where rt describes 
the transition radius between the rising and flattening of 
the curve), or a multi-parameter fit with rt = 2.1kpc, 
7 = 1.5. To compare the shape of the rotation curve 
with Courteau's sample we compute ri/ropt (where ropt 
is the radius enclosing 83% of the light in the source plane 
photometry) to be ~ 0.25. This ratio is higher than most 
galaxies in the local sample, but not anomalous. Around 
10% of local bright galaxies have similarly shaped rotation 
curves. 



^ Based on observations made with the William Herschel Telescope operated in the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. 
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Fig. 3. — Examples of the [On] emission doublet as seen in the 
reduced IFU data. Each line of the image is the spectrum from a 
single IFU lenslet. These are arranged in blocks, each block cor- 
respond to an East- West slice through the short dimension of the 
galaxy. The rotation of the galaxy can be seen in each block as well 
as between blocks, clearly showing that the lensed galaxy's dynam- 
ics are resolved in both dimensions. 

3. DISCUSSION 

The tight correlation between luminosity and rotation 
velocity for spiral galaxies in the local universe is known 
as the TuUy-Fishcr (TF) relation (TuUy & Fisher 1977). 
With a single high redshift galaxy, we can measure the 
evolution of the offset in the TF relation under the as- 
sumption that the slope remains fixed. In Fig. 4, we com- 
pare the rotation velocity and source brightness of arc # 
289 with that of local and other high redshift galaxies in 
rest frame B and I. Whilst the current high redshift data 
are concentrated in rest frame B-band, rest frame 7-band 
observations will prove a more rigorous test of evolution of 
the TF relation since the corrections for dust and on-going 
star formation are much smaller at longer wavelengths. 
The rest frame /-band TF therefore gives a clearer indica- 
tion of the true stellar luminosity and hence the ratio of 
stellar mass to total halo mass. 

The position of the galaxy on the TF relation in both 
the B and / bands shows good agreement with local data 
(Pierce & Tully 1992; Mathewson & Ford 1992; Haynes et 
al. 1999; Tully & Pierce 2000; Verheijen 2001). The galaxy 
is slightly offset to brighter magnitudes in the _B-band re- 
lation, as we would expect from its blue colors. This data 
is in agreement with existing intermediate (z ~ 0.5, Vogt 
ct al. 1997; Milvang- Jensen ct al. 2001; Bohm et al. 2002; 
Zieglcr ct al. 2003) and high redshift studies {z ^ 1, Vogt 
et al. 2002; Barden et al. 2003). It should be noted, how- 
ever, that these studies generally have to infer the turn 
over in the rotation curve from deconvolution of the rising 
rotation speed and the seeing. This might lead to a sys- 



tematic underestimate of the asymptotic rotation speed. 
In contrast, our data clearly resolve the shape of the ro- 
tation curve (Fig. 4 inset) leading to much smaller uncer- 
tainty in the rotation speed of an individual galaxy. 

Clearly, it would be dangerous to draw far reaching con- 
clusions from a single galaxy, but the following comparison 
illustrates the sensitivity that can be expected if this level 
of evolution is confirmed by further observations. The the- 
oretical evolution of the J-band TF relation from hierarchi- 
cal models of galaxy formation from Cole et al. (2000) prc>- 
dict that for any given disk circular velocity, the / band lu- 
minosity should decrease by ~ 0.1 magnitudes from z = 
to z = 1, whilst in the B-band such models predict an 
increase in luminosity of ~ 0.5 magnitude for the same 
redshift change. In contrast, we would expect a decrease 
in luminosity of ^ 0.7 magnitudes from z = to 2 = 1 if 
all of the galaxy's mass were already in place at z = 1, but 
only half of the stars had yet formed. The small change in 
the /-band magnitude shown by arc#289 suggests a prefer- 
ence for hierarchical rather than the "classical" formation 
model. 

The methods we have developed in this paper show the 
potential of IFU observations of gravitational lensed galax- 
ies as a means of studying the spatially resolved properties 
of high redshift galaxies in a remarkable level of detail. 
While we have concentrated here on the dynamics of the 
galaxy #289 in A 2218, our data can also be used to in- 
vestigate the; spatial distribution of star formation. By 
combining these optical data with observations with near- 
infrared integral field units, it will be possible to study 
the emission line ratios of [Oii]:[Oiii] and Hq::H/3. We 
could then determine the distribution of reddening across 
the galaxy and its spatially resolved chemical abundance. 
This would provide powerful insight into the nature of star 
forming galaxies at this early epoch. 
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Fig. 4. — Left: ylrc#289 on the Tully-Fisher relation in rest frame S-band compared to high redshift {z ~ 0.83) field galaxies 
(Milvang- Jensen et al. 2001) and the high redshift sample from Vogt et al. (1999). For comparison we show the low redshift local 
fit from Pierce & TuUy (1992). The solid triangle shows a massive disk galaxy (L451) at 2; = 1.34 (van Dokkum & Stanford 2001). 
The galaxy rotation curve (inset) shows the peak-to-peak rotation velocity of the arc in [On] emission built from the IFU imago 
in the source plane. The error bars shown are formally 3-a and alternate points are independent. The horizontal error bars show 
0.7" seeing transformed to the source plane. Right: The rest frame /-band Tully-Fisher relation compiled from Mathewson & 
Ford (1992) and from the Ursa-Major Cluster (Verheijen 2001). Arc#289 is shown by the solid point and lies very close to the 
mean TF relation for present day spirals. The small change in the /-band magnitude shown by arc#289 suggests a preference for 
hierarchical rather than the "classical" formation model. 
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